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3^ ■ Abstract 

0\ . 

■ There are now several analyses reporting quantized differences in the 

OO ' redshifts between pairs of galaxies. In the simplest cases, these differential 

redshifts are found to be harmonics of fundamental periods of approxi- 
1 mately 72 kms -1 and 37.5 kms -1 . In this paper a wave equation is 

derived based on cosmological general relativity, which is a space- velocity 
O ■ theory of the expanding Universe. The wave equation is approximated 

to first order and comparisons are made between the quantized solutions 
5^ ' and the reported observations. 

99' 

; 1 Introduction 

^ ' There have been numerous, detailed observations and analyses of discrete dif- 

ferential red shift found between pairs of galaxies, especially TifFt PQ| and Napier 
& Guthrie|2]. Attempts have been made in the past by Cocke[3| and Nieto|H 
to develop a gravitational wave equation which could explain the results. 

In this paper a linear second-order wave equation is derived based on the 
space- velocity theory of cosmological general relativity (CGR) of Carmeli[S]. A 
spherically symmetric metric is used with comoving coordinates. First order 
approximations are made for small and large radial distances and the resulting 
equations are solved. The large distance quantized solutions are compared to 
the observational data from differential redshifts found between galaxies. 
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2 Wave Equation 



The CGR spherically symmetric, comoving space- velocity metric (Ref. Eq. 
(A. 5)) is defined by the line element 

ds 2 = r 2 dv 2 - e^dr 2 - Q 2 (d0 2 + sin 2 {6)d<f) 2 ) , (1) 

where coordinate v is the radial velocity of expansion of the universe, and r 
is the Hubble-Carmeli time constant, its value is r = 12.486 Gyr (Ref. [5], 
Eq. (A. 66), p. 138). The functions p and Q are dependent only on the radial 
velocity v and the radial coordinate r, where r, 9 and <f> are the usual spatial 
spherical coordinates. From Eq. ^the non-zero elements of the metric g^ are 

5oo = 1, (2) 

.911 = -e", (3) 

.922 = -Q 2 , (4) 

533 = -Q 2 sin 2 {6). (5) 

A postulate of this cosmological theory is that the metric g^ v satisfies the 
Einstein field equations 

Gfj,u = Rfiu — -^9fj,uR = ' (^) 

where i? M „ is the Ricci tensor, R = g a ^R a p is the scalar curvature, Tn V is 
the momentum-energy tensor, n = 8itG/(c 2 t 2 ) , G is Newton's constant and 
c is the speed of light in vacuo. (Compare this to the Einstein equation in 
General Relativity theory where k — 87rG/c 4 .) The momentum-energy tensor 
of a perfect fluid (Ref. 0, Eq. (A.10)) is 

Tp U = Peff u^Uy + p (u^u u - g^) , (7) 

where the effective mass density p e ff = p — p c , where p is the average mass 
density of the Universe and p c is the critical mass density, a constant in CGR 
given by p c — 3/(87tGt 2 ). Also, p is the pressure, and u M is the four- velocity 

u a = u a = (1,0,0,0). (8) 

The solution for Eq. E|was derived (Ref. @, Appendix A. 4) with the results, 

Q = r, (9) 

= eMr) = TTTW' (10) 

f(r) = l ~J^r 2 , (11) 
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where the density parameter £1 = p/p c - 
The linearized line element is defined by 

ds = d (tv) — ie^'^dr — irdO — ir sin(9)d(p , 



(13) 



where i = 

The linearized first order spinor differential equation follows by canonical 
correspondence [7] from Eq. 1131 using spinors , 
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Square Eq. 1 141 and assume the r.h.s. linearizes to the form of the d'Alembertian, 
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In CGR the condition for the expansion of the universe is defined by setting 
ds = 0. For the wave equation, it is assumed that the expansion of the universe 
corresponds to setting d 2 ^ /ds 2 — 0. With this condition, Eq. 1 151 becomes 



t 2 dv 2 r 2 dr 
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(16) 



3 Solving the Wave Equation 

We solve Eq. by a separation of variables. Assume that fy(v,r,9,<fi) = 
^o(v)^i(r)^ 2 (S)'^3(4>)- Substituting the composite function for "J" into Eq. 
1161 and dividing both sides by iff we obtain 

1 d 2 ^ (v) _ 



if>o(v) r 2 dv 2 

1 1 d 



-M(r)_ 



if>i(r) r 2 dr 



1 
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d9 
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(17) 
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In Ea ll7l since the left hand side is a funtion of (v) only, while the right hand 
side is a function of (r, 9, <fi) only, they must each equal a constant, — D 2 . Thus 
we can say 



1 d^ (v) = _ d2 



D 2 = e-"W 



1 



1 Id 



(18) 



\l/i(r) r 2 dr 
1 d 



V 2 (9) sin{6) dQ 

1 1 



sin(9) 



irfWi(r) 
dr 

d^ 2 (9) 
dd 



d 2 ^ 3 (</>) 

Solutions for Eq. have the form 



(19) 



(20) 



where D is called the intrinsic curvature. 

Eauation llijl can be split in two parts, one a function of (9, <fi) only, the other 
a function of (r) only, both equal the same constant — L 2 . They simplify to 

1 Id/. « 

L 2 = 



and 



* 2 (#) sm(0) d0 

1 1 d 2 fr 3 (</>) 
' sin 2 (9)^ 3 (<f>) dej) 2 



(21) 



J_ d_ 

r 2 dr 



dr 



+ £>V (r) *i(r) ^— *i(r) = 0. 

The solutions for Eq. 1211 are the well known spherical harmonics 
* 2 (0)* 3 fo) = -PftccM^e*** , 
where Pi(cos(9)) are the associated Legendre functions and 

k = 0,±1,±2,±3,..., 

L 2 = 1(1 + 1) , 



(22) 

(23) 

(24) 
(25) 



l=\k\, \k\+l, \k\+2, 



(26) 
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To solve Eq. [221 for 'l'i(r) we make first order approximations for e^ r ^ 
Assuming | (1 — Q) r 2 / (c 2 t 2 ) 1 we have from Eqs. ITU1 and ITT1 



(l-fi)r 2 



(27) 



Putting this approximation into Eq. 1221 and collecting terms we have 



f d ( 2 d*i(r)\ L 



r 2 



dr 
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2 , ( 1 ~") I 2' 



*i(r) = 0. 



(28) 



4 Approximation for Small r, Quantization 

For small distance r from a comoving source, such as a galaxy center, and for 
small L 2 , in Eq. |^H1 drop the third term and drop the second of the forth term, 
and subtitute L 2 = 1(1 + 1) to obtain 



1 d 

r 2 dr 



id$i(r) 
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2 ^+1) 



*i(r) = 0. 



(29) 



Eauationl29lis identical in form to that of Ref. 0], (Eq. 5.2). The eigenfunctions 
are the two linearly independent spherical Bessel functions ji(r D) and yi(r D), 
with eigenvalues determined at the boundary r = cr where ji(crD n i) = and 
yi(crD nl ) = for n = 1,2,3,4,... and I = 0,1,2,3,.... The eigenfunctions are 
related to the ordinary Bessel functions Ji(x) by 



3l{x) = \/ 2^+1/2(2;) 



(30) 



wW = (-1) 



J_i_i/2(a;) • 



(31) 



The first 6x6 values of the zeros crD n i of jt (x) are 
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1 


3.14 


4.49 


5.76 


6.99 


8.18 


9.36 


2 


6.28 


7.73 


9.10 


10.41 


11.70 


12.97 


3 


9.42 


10.90 


12.32 


13.70 


15.04 


16.35 


4 


12.57 


14.07 


15.51 


16.92 


18.30 


19.65 


5 


15.71 


17.22 


18.69 


20.12 


21.53 


22.90 


6 


18.85 


20.37 


21.85 


23.30 


24.73 


26.70 



(32) 
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For high order n > N, the zeros of ji(cTD n i) approach 

lim / l\ 

crD rd = [n + - pr , (33) 



n — > oo \ 2 

n = AT,iV + l,iV + 2,iV + 3,..., (34) 
I = 0,1,2,3,.... (35) 

The discrete intrinsic curvature D n i are due to the gravitational wave modes 
associated with the expansion of the Universe. It is assumed that the intrinsic 
rcdshift Zi observed in the wavelength of a photon is proportional to the intrinsic 
curvature of the source of the photon, 

z i<n i = bcrD n i, (36) 

where b is a dimensionless positive constant. These eigen solutions are for 
small distances from some comoving source. Notice that under the prevailing 
assumptions these solutions are independent of f2, which means that they apply 
to all three phases of the expansion of the Universe: deceleration, constant 
and accelerations. A prediction of CGR is that the expansion is currently 
accelerating Pj . 

For / = or for large n, Eq. OHHdoes yield a linear relationship of intrinsic 
redshift to intrinsic curvature. However, these are for small distances, whereas 
in the reports to be discussed later on, the galaxies can be separated by large 
distances of a half megaparsec to a few tens of megaparsec. We suggest that it 
may be applicable in the analysis of the redshift of quasars presumed to be in 
near proximity to a host galaxy as reported in Ref. HD and |12| . 



5 Approximation for Large r 

For large r make the approximation 



1 d 

r 2 dr 



d^ 1 (r)\ d 2 *i(r) 



dr 



dr 2 



For large distance r from a comoving source and for small L 2 , in Eq. 
the second term and the second of the fourth term to obtain 



d 2 ^i(r) 
dr 2 



D 2 _ (LJ?) D 2 r 2 



*i(r) = 0. 



(37) 
ignore 

(38) 



There are three possible solutions to consider, depending on whether fi > 1, 
n = lorO<fi<l. 
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5.1 > 1, Decelerated Expansion 

For > 1, generalized power series solutions are required for Eq. |2H1 A solution 
for ^i(u) has the form 

oo 

*i(u) - ^afcA (39) 

fe=0 

u = Voir, (40) 



where a = VSl — ID/ (cr). When ^\{u) is substituted into Eq. OH the values 
for the coefficients are 

oa = (41) 
as = (42) 

ak+2 ~ (fc + 2)(* + l) ' (43) 

fc - 2,3,4,..., (44) 

where 7 = Dcr/V^ — 1. We can let ao = 1 and a\ = 0. Thus all the odd 
numbered coefficients a,2k+i vanish. The condition for the series Eq. to 
converge absolutely by the ratio test is 

lim|(^±^)|=p, (45) 
fe^oo \ a 2 k J 

where p < 1. There is apparently no need for any discrete values of the param- 
eter D. 

5.2 f2 = 1, Constant Expansion 

For — 1 , Eq. [38] has the solutions 

= 9 c e kiDr . (46) 
These solutions allow continuous values for D. 

5.3 < Q < 1, Accelerated Expansion, Quantization 

For < ft < 1, Eq. 1381 is the equation of the quantum harmonic oscillator with 
solutions involving the Hermite polynomials H n {\) of order n, 

*i(x)=*c(V x2/2 )#n(x), (47) 
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where \ — (1 — SI) 1 / 4 (D/ct) 1 ^ 2 r, and ty c is a constant. The solutions have the 
requirement that 

ctD„ 

^== = ^+1, (48) 

where n = 0, 1, 2, 3, .... 

As stated above, the intrinsic redshift Zi is assumed to be proportional to 
the intrinsic curvature D ni 

Zitn = bcTD n = Z (2n+l) , (49) 



z Q = bVT^n, (50) 

where 6 is a dimensionless positive constant. 
5.4 Theory Versus Experiment 

We will now discuss the quantized version of the solutions for large r given by 
Eq. This equation is valid during the accelerated expansion stage. The idea 
put forth here is that the wave equation governed the accelerated expansion 
of the Universe and created states of quantized intrinsic curvature in the fluid. 
When the fluid coalesced to form gas clouds, the quantized curvature states 
persisted. 

Assume that 6=1. From Eq. [5U|we get 

z$ = i-n. (51) 

A combined redshift formulajlO) will show the relationship of the intrinsic red- 
shift to the total redshift z t , 

1 + Z t = (1 + Z C ) (1 + Z r ) (1 + Zi) , (52) 

where z c is the cosmological redshift due to the Hubble expansion and z r is a 
random component due to all other effects. From Eq. [521 the difference of total 
redshift between a galaxy in intrinsic state ri2 and another galaxy in intrinsic 
state m, where n\ < «2, is 

Zt,2 - Z t ,l = (1 + Z c . 2 ) (1 + Z r fi) (1 + Zi,n 2 ) 



-(l + z c , 1 )(l + z r>1 )(l + z i>ni ) . (53) 

Assume that the cosmological components z c , m are small and approximately 
equal for both galaxies and that their random components are negligible. Com- 
bining Eqs. 0!51and|nnia,nd retaining only the first order terms we get differential 
redshift 

A.Z fc = z t ,2 - z t ,i « Zi,„ 2 - z i>ni = 2 Z k , (54) 
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where k = n% — n\. From Eq. 1541 define the differential intrinsic velocity 

AV k = cAZ k = 2cZ k. (55) 



Now, we will briefly discuss two reports on the quantization observed in 
the differential redshifts from normal galaxies. In Ref. pQ, 31 pairs of galaxies 
were observed in radio and some in optical. Their redshift differentials were 
found to have a fundamental periodic velocity of AV = 72kms _1 . We use 
c = 299792.456 kms" 1 for the vacuum speed of light. To reproduce these peri- 
odicities, set Zo = 36/c = 1.201 x 10 -4 . From Eq. this would correspond to 
a density parameter 1 — CI = 1.442 x 10 -8 . From Eq. [55]we have 

AV k = 2cZok = 72kkms-\ (56) 
k = 1,2,3,4. (57) 

Equation EfH gives the sequence of periodic velocities 72, 144, 216 and 288 
[kms -1 ] that appear as peaks in (Fig. 1) of Ref. P^. 

In Ref. [5], 97 galaxies were observed, where the velocity of the galaxies 
extended out to 2600 km s -1 . Their redshift differentials were determined, rel- 
ative to the Galactic center, with a fundamental periodic velocity of AV = 
37.5 kms" 1 obtained with amazing regularity. Refer to their (Fig. 4) of that 
report. To reproduce these periodicities, set Zq = 18.75/c = 6.2543 x 10 -5 . 
This would correspond to a density parameter 1 — fl = 3.9116 x 10~ 9 . From 
Eq. [55] we have 

AV k = 2cZ a k = 37.5 k kms -1 , (58) 
k = 1,2,3, ...,38. (59) 

It is evident that Eq. [SSI gives the sequence of periodic velocities that appear 
as peaks in (Fig. 4) of Ref. 0, from 37.5 kms -1 up to 1425 kms -1 . 

In closing we note that these are only preliminary analyses to show the 
efficacy of quantized intrinsic redshift. A more thorough analysis would study 
galaxy populations in terms of cosmological redshift and parametric density. 

The author is grateful to the Joint Astronomy Centre, Hilo, Hawai'i for 
many years of employment and support, and to Professor Moshe Carmeli for 
his theories and communications. 
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